Microwave conductivity of a d-wave superconductor disordered by extended 
impurities: a real-space renormalization group approach 
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Using a real-space renormalization group (RSRG) technique, the microwave conductivity of a 
d-wave superconductor disordered by extended impurities is calculated. To do this, a semiclassical 
approximation is invoked which naturally accesses the Andreev bound states localized near each 
impurity. Tunneling corrections (which are captured using the RSRG) lead to a derealization of 
these quasiparticles and an associated contribution to the microwave conductivity. 
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I. INTRODUCTION 

One aspect of the electronic properties of the 
cuprate superconductors which remains mysterious is 
their frequency-dependent microwave conductivitji^, 
the measurement of which sheds light on the effect of 
disorder on quasiparticle properties. A common fea- 
ture of many theoretical approaches to disorder (see, 
e.g., Refs.y0) m d-wave superconductors (although there 
have been a number of exceptions£*§*£*&) is the assump- 
tion that the scattering may be modelled by pointlike 
impurities. Recent work 7 found that the spectrum of 
quasiparticle excitations is qualitatively different when 
the disorder is taken to consist of extended impurities, 
i.e., impurities characterized by a typical size a Af, 
where Af is the Fermi wavelength. While the density 
of states p(E) is believed to vanish linearly at low ener- 
gies for the case of poinlike impurities&iS, in Ref.0it was 
found that, for disorder consisting of extended impurities, 
p{E) is divergent (~ l/_Elog 3 i5) at low-energies. This 
low-energy buildup of states arises from the hybridiza- 
tion of Andreev bound states occurring near each such 
extended impurity. 

Before proceeding, we provide two specific experimen- 
tal motivations for studying the effect of extended im- 
purities on the electronic properties of d-wave supercon- 
ductors: Firstly, recent technical advances have allowed 
the fabrication of ultra-pure YBCO which has almost no 
atomic disorder in the Cu02 planes 2 . A significant source 
of disorder may originate away from the Cu02 planes and 
can be expected to project a relatively long-wavelength 
potential on the CuC>2 planes. Secondly, there is also 
the possibility of deliberately creating extended defects 
via ion irradiation techniques 11 . Apart from these spe- 
cific possibilities, we note that although the calculation 
presented here strictly applies to the case of disorder sat- 
isfying a > Ap, the physical effects may persist beyond 
this regime. 

The purpose of this Paper is to further explore the 
consequences of extended impurities on the electronic 
properties of a d-wave superconductor by computing the 
associated contribution to the microwave conductivity. 



To accomplish this, we invoke a semiclassical approxima- 
tion&L^ which accesses the low-energy Andreev-bound 
states occuring near each impurity. Although the indi- 
vidual Andreev bound state wavefunctions are localized 
near the impurities, they become delocalized at low ener- 
gies due to tunneling corrections^, leading to a concomi- 
tant contribution to the conductivity. We find that a 
natural way to study such tunneling effects between An- 
dreev bound states (and to compute their associated con- 
tribution to the conductivity) is via a real-space renor- 
malization group (RSRG) technique which is a variant 
of the decimation technique familiar from random spin 



systems' 
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Thus, although the Andreev bound 



states near individual impurities are strongly coupled to 
each other, the RSRG reveals effective Andreev bound 
states associated with many impurities which are weakly 
coupled and which contribute strongly to the low-energy 
microwave response. 

This paper is organized as follows. In Sec. ^ we em- 
ploy a semiclassical approximation which reduces the 
Bogoliubov-de Gennes eigenproblem for a disordered d- 
wave superconductor to a family of one dimensional ran- 
dom pair potential models (represented by a Hamiltonian 
Ha). In Sec. IIHI we introduce the RSRG decimation 
procedure for computing average properties of Ha- In 
Sec. IIVI we apply this procedure to the computation of 
the microwave conductivity of a d-wave superconductor 
and discuss our results. In Sec. E] we make some brief 
concluding remarks. In the Appendix, we present a cal- 
culation of the average single-particle Green function of 



II. QUASIPARTICLES OF A DISORDERED 
D-WAVE SUPERCONDUCTOR 



We assume that the quasiparticles of a d-wave super- 
conductor with extended impurities are governed by the 
following Bogoliubov-de Gennes (BdG) action: 



S = 



dr / d 2 r¥ 



h A 
A -h 



(1) 



2 



where h = — V 2 — e F + V(r), ep = Pp is the Fermi energy 
[i.e. pf (= 27r/Ap) is the Fermi wavevector] , and we have 
adopted units in which h 2 /2m = 1. The d-wave pair 
potential operator A = A (p 2 — Vy)lv\ with A being the 
pair-potential maximum. We take the disorder potential 
V(r) to arise from randomly located impurities of typical 
size a and typical spacing £, and focus on the regime 
i > a. 

To study the Andreev bound states occuring near each 
extended impurity, we next invoke a semiclassical approx- 
imation (valid for pfo > 1) which has been extensively 
discussed elsewhere&LiS and which exchanges the BdG 
Hamiltonian for an effective problem residing on a clas- 
sical trajectory r c (s) = (x c (s),y c {s)) solving Newton's 
equation [i.e. 2p|9 2 r c (s) = — VV(r c (s))]. This is imple- 
mented by writing for a particular incoming momentum 
direction n (— r c (s) for s — > -co) the field \& as 

* = Ae ipFS ^, (2) 

and expanding to leading order in derivatives. Here, the 
eikonal S satisfies VS* = r c (s) with the overdot denoting 
differentiation with respect to the parameter s along a 
trajectory. The amplitude A ~ 1. A given realization of 
V(r) leads to a set of classical trajectories labelled by n 
and an impact parameter b. 

By inserting Eq. (0) into Eq. JQ), we see that the dy- 
namics of the trajectory-dependent quasiparticles (rep- 
resented by the field ip) is governed by the Andreev 18 
Hamiltonian Ha 

_ (-2i PF d s A(s) \ 
Ma ~\ A( S ) 2ip F d s )> W 

which has the form of a one-dimensional superconduc- 
tor with an effective pair potential A(s) = Ao(a; c (s) 2 — 
y c (s) 2 ). For a given classical trajectory, A(s) exhibits 
rapid variations that we wish to treat via an averaging 
procedure. To properly compute the disorder-averaged 
correlators of A(s) one must solve Newton's equation for 
r c (s) in the presence of a given realization of V(r) and 
then average over all such realizations. We shall not at- 
tempt such a difficult calculation; however, on general 
grounds we expect that A(s) has zero mean and short- 
range correlations beyond £. We shall analyze Ha by 
studying the Andreev bound states associated with sign 
changes in A(s) using the RSRG. Before discussing how 
this works, let us briefly review previous results for this 
one-dimensional random pair-potential Hamiltonian (we 
note that Ha has appeared in several other condensed- 
matter contexts 19 ). The disorder-averaged low-energy 
density of states p(E) ~ 1/E\v? E of Ha was first ob- 
tained by Ovchinnikov and Frikhman 20 using a Fokker- 
Planck equation. In Ref. |3, this same result was ob- 
tained by mapping Ha to a one-dimensional random hop- 
ping model. More recently, the disorder-averaged single- 
particle Green function G(x — y; E) of Ha was obtained 
by Balents and Fisher— using a supersymmetry method. 
They found that G(x — y; E) ~ exp[— \x — y\j log 2 E] for 



| a; — y\ — > oo. In the Appendix, we reproduce these re- 
sults using the RSRG technique; here we note that, in 
the present context, this behavior signifies that the An- 
dreev bound state wavefunctions are exponentially local- 
ized but with a characteristic length scale that diverges 
at low energies (i.e. they become delocalized). Finally, as 
discussed in Refs00, Ha is closely related to models of 
supersymmetric quantum mechanics 22,23 . 



III. REAL-SPACE RENORMALIZATION 
GROUP APPROACH 

In this section we discuss how the disorder-averaged 
low-energy properties of Ha may be analyzed via a RSRG 
approach. Of central importance in such an analysis are 
the sign-changes (i.e., zeroes) of A(s) 24 . Henceforth, our 
qualitative picture of A(s) along a trajectory is depicted 
in Fig.^ i.e., it exhibits many sign changes. Before con- 
sidering the many sign-change case, however, let us con- 
sider an infinite system in which A has exactly one sign 
change. In this case, Ha has an exact zero-energy (ZE) 
eigenstate of the ioiu&^^ 

^WKf^e^^^W, (4) 

where the +(— ) corresponds to the case A(s) < 
(A(s) > 0) for s > s n and we have omitted an overall nor- 
malization factor. Indeed, the normalizability of Eq. Q 
relies on the existence of the sign change in A(s). The 
subscript n enters when we consider the case in which 
the pair potential has multiple sign changes; henceforth 
it shall refer to the nth zero s n of A(s) and refer to ^ n (s) 
as the ZE state associated with s n . In the case in which 
A(s) has multiple sign changes, however, we must ac- 
count for tunneling corrections t n = (^ n \HA\^ n +i) be- 
tween nearest-neighbor ZE states (i.e. they are no longer 
at zero energy). A direct calculation gives 6 

*» = \l^\&'(Sn)&'(s n+1 )\ 1 / i (5) 

xexp (^ir +iA(s)ds i)' 

where the prime indicates differentiation with respect to 
s. 

Having identified the ZE states ^ n associated with the 
zeroes of A(s), we now turn to the RSRG analysis of 
such states. Although it is formally similar to that stud- 
ied in Refs. there are conceptual differ- 
ences which are important. The essential result of such 
an analysis is the following: Although the matrix ele- 
ments t n coupling the may be large, there are effec- 
tive ZE states associated with many sign changes which 
are weakly coupled and which may be treated via pertur- 
bation theory. Furthermore, the couplings between such 
effective ZE states have a probability distribution which 
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FIG. 1: Sketch of the effective pair potential along a trajec- 
tory. The pair-potential has been deliberately drawn smaller 
between sites S2 and S3 to illustrate a point; see the text for 
details. 

is universal in a sense discussed below. To begin, we note 
that for a given realization of the pair potential, there ex- 
ists a strongest matrix element i max . The ZE states as- 
sociated with this matrix element are strongly coupled, 
forming symmetric and antisymmetric wavefunctions at 
energies ±f max . The difference between these energies 
defines our initial bandwidth Qq = 2i max . To estimate 
the magnitude of Oo, we first assume the impurities are 
dense enough such that the exponential factor in Eq. (0 
is negligible. Then, we assume that A'(s„) « Ao/a near 
an extended impurity of size a. Using Ao ~ 50meV, we 
have f2o ~ 10 14 s~ 1 / v / pf«- 

Generally, we are interested in an observable at some 
low-energy lu <C fio; the RSRG progressively elimi- 
nates all pairs of ZE states coupled by matrix elements 
larger than u>. To see what this means, let us examine 
Fig. H According to Eq. 10, the matrix elements t n are 
exponentially dependent on the pair-potential integral 
I f° n+1 A(s)as| between pairs of sites. Thus, in Fig. 2] the 
strongest matrix element is that connecting sites S2 and 
S3, (i.e., t%) since the aforementioned integral is smallest 
for them. To obtain the decimation procedure, we recall 
that the the sites s n are defined by the sign changes of 
A(s) along a particular trajectory. The fact that A is 
only slightly above zero for the segment between S2 and 
S3 suggests that we could have simply ignored the sign 
changes at S2 and S3 without making any appreciable er- 
ror in computing physical quantities. This defines a new 
description of Ha- a basis of ZE states that are associ- 
ated with all sign changes of A(s) except those at sites 
S2 and S3. The effect of these zeroes of A(s) is incorpo- 
rated into a weak matrix element between the effective 
ZE states at sites si and S4 (that are each associated with 
three sign changes of A(s)). This procedure exchanges 
four ZE states for two weakly-coupled ZE states (and 
lowers the bandwidth to £1 < flo). The associated ma- 
trix element t between the effective ZE states at sites si 
and S4 is still given by the general formula Eq. JSJ; it is 
easy to see that this formula has the following recursion 
property: 

t = ht 3 /t 2 . (6) 

Remarkably, this recursion relation is the same as that 
obtained by Fisher— for the random antiferromagnetic 



spin chain; the existence of such a relation is the reason 
the RSRG works in this seemingly unrelated system, ft 
is convenient to define (as was done in the random spin 
chain casei^) logarithmic matrix-element strengths ( n = 
log(£l/2) — logt„ and arc lengths l n = s n+ i — s n . It is 
useful to note that up to a constant (defined so that the 
minimum is zero), £„ is essentially the abovementioned 
pair-potential integral | J Sn+1 A(s)ds| that motivated the 
decimation scheme. A particular (trajectory-dependent) 
form for Ha defines an initial probability distribution 
P(C, I', ^0) for the {Cn} and {/„}, in which ZE states have 
been assigned to all sign changes of A. How do these 
values evolve under the RSRG decimation? Once sites 

52 and S3 have been decimated, the new effective arc 
length is clearly I = S4 — si. This, along with Eq. JfjJ), 
implies that 

C = Ci+C3-C2 = Ci+C3, (7a) 
I = h+h + k, (7b) 

where in Eq. (|7a|l we have used the fact that (as noted 
in Ref. Il4f) since the matrix element connecting S2 and 

53 is assumed to be the strongest, £2 = by definition. 
By iterating this procedure, we produce a basis set of 
ZE states which are more weakly coupled each iteration 
(and which are associated with progressively more sign 
changes of A(s)) but which still contain all the low en- 
ergy physics of Ha- As mentioned above, although their 
physical origin is slightly different, the recursion rules 
for the {£„} and the {/„} are the same as those studied 
by Fisher for random spin chains. The essential results 
of this work (for our purposes) are the number of sites 
710 = 1/ In 2 Qq remaining after having lowered the 
bandwidth (by decimating) to energy f2 and the asymp- 
totic (J7 — > 0) probability distribution P((,l;fl), which 
is given byi^ 

P((,i-,n) = ^Q(c/r,//r 2 ), (8a) 

Q(r],y) = _AA^_ e -"^ coth ^. (8b) 
sinh v cy 

Here, Q(r], y) is the Laplace transform of Q(rj, y) and T = 
ln$7o/^- The integration constant c defines the length 
scale over which I varies. We assume this is given by the 
typical value of / in the initial distribution, i.e., c ~ £. 

IV. MICROWAVE CONDUCTIVITY 

Having outlined the RSRG technique, we now apply it 
to the calcuation of a specific quantity: The microwave 
conductivity of a d-wave superconductor disordered by 
extended impurities, averaged over all realizations of the 
disorder. Before embarking on the calculation, let us 
discuss the physical picture. We are interested in the 
possibility of the Andreev bound state quasiparticles con- 
tributing to the microwave conductivity. From the outset 



4 



it is not obvious that they will do so, since they are nomi- 
nally localized [see Eq. Q] near a particular sign change. 
The results of Balents and Fisher for the Green function 
G(x — y; E) of Ha give hope to this possibility revealing 
that the quasiparticles become delocalized at asymptoti- 
cally low energies. As the conductivity is related to a two- 
particle Green function (see below), the results of Ref . l2ll 
are not sufficient to compute the conductivity. However, 
as we shall see the peculiar way in which tunneling cor- 
rections conspire to form effective Andreev bound states 
associated with many sign changes conspires to give an 
interesting contribution to the microwave conductivity. 

For a given realization of the disorder, the real part 
of the microwave conductivity o\ is given by the Kubo 
formula 25 : 



<Ji{lo) = --Imn R H, 



(9a) 



Il(ifl) = -^ /3 dre^<j(q,T).j(q,0))| q ^ o , (9b) 

Here, A is the area of the system and to obtain the re- 
tarded polarization tensor we use the usual analytic con- 
tinuation n R (ti>) = H(ifl — > u> + iO + ). The angle brackets 
represent the thermodynamic average with respect to the 
action S. We shall denote the disorder average by square 
brackets; however, before attempting to evaluate this av- 
erage we first simplify Eq. for n(iO). 

Within the semiclassical approximation of Sec. [HI 
j(q, t) is obtained by inserting Eq. into the usual 
current operator j(r, r) = ie(^^V^ — /i.e.). The Fourier 
transform j (q, r) is (in the q — ► limit) 

j(0,r) = -epll] / dbdsr c (s)^(s,T)TP(s,T). (10) 



Standard manipulations then yield the following expres- 
sion for (Ji (lj): 

cti(co) = — / dbdsds'v c {s) ■ Y c (s')lmK K (s,s';uj), 

n 

(11a) 

K(s,s';t) = (^(s,t)iP(s,t)^(s',0)iP(s',0)), (lib) 

Next, we use Wick's theorem (S is a Gaussian action) 
to express the Matsubara Fourier transform K(s,s';i£l) 
of Eq. (|llbjl as a sum of two terms, one of which gives 
a vanishing contribution to H The other contribution 
to K(s, s'; in) may conveniently be expressed in terms 
of the normalized eigenstates {ipi(s)} (with eigenvalues 
{£.,}) of ff A : 

x(n F (Ei) - n F (Ej))S(uj + E t - Ej), (12) 

where, since we are interested in ImII R (w), we have dis- 
played lmK R (s, s'; u>) directly. 



Turning to Eq. Ijllafl . we see that for the conductivity 
we require the disorder average of lmK R (s, s'; ui) mul- 
tiplied by r c (s) • r c (s'). At low energies, the ZE states 
hybridize and become delocalized along the classical tra- 
jectories. By contrast, the average of r c (s) ■ i" c (s') is ex- 
pected to vanish on long length scales as the classical 
trajectories bounce off the extended impurities. Thus, 
it is sensible to approximate the average of the product 
of these quantities by the product of the averages; fur- 
thermore we shall approximate for the classical problem 
[r c (s) -r c (s')]dis. — exp(— |s— s'\/t) i.e., exponentially de- 
caying correlations. We believe that the precise form of 
this correlator is unimportant at low energies. However, 
we note that the virtue of this form is that it naturally in- 
terpolates between ballistic motion on short length scales 
(i.e. [(r c (s) — r c (0)) 2 ]dis. = s 2 for s — > 0) and diffusive mo- 



Tc(0))< 



2s/e 



tion on long length scales (i.e. [(r c (s 
for s — ► 00). 

Next, we turn to the computation of the disorder av- 
erage of lmK R (s, s'; u>), for which we require the RSRG 
analysis. The purpose of the RSRG is to eliminate (dec- 
imating until Q = lo) all the strong matrix elements so 
that we are left with an effective theory having only weak 
(i.e. t n < lo) matrix elements between ZE states, allow- 
ing the use of perturbation theory. Following Fishei~, 
we assume that the broadness of P(C, l] fi) indictates that 
the disorder average is dominated by a class of rare pair- 
potential configurations that provide a large contribu- 
tion. For the disorder average of Eq. 112(1. the class we 
have in mind is when an undecimated pair of states is 
separated by Z = \s — s'\ exactly; the relative number 
of such configurations with logarithmic matrix element 
strength £ is simply given by n w P(C, \ s — s'\;u>). For a 
given £ and to leading-order in perturbation theory (a 
similar calculation was done in Ref. ITtT) . the ZE states 
at s and s' form symmetric (S) and antisymmetric (A) 
wavefunctions with energies ±t = ±^e - ^ (recall we have 
set Q = of). At low energies, the only terms which con- 
tribute appreciably to the sum in Eq. 1|12[> are these states 
(i.e. i,j = S,A). Thus, we have (upon summing over all 
such configurations which amounts to integrating over £) 



[lmK R (s,s';uo)} dis . cx n w / d£P{C, \s - s'|; u) (13a) 

Jo 

x6(oj - we _c )tanh(w/4T), 

aP(0,|— '|;o,)^^3b) 
a; log llo/u) 

where we have omitted overall temperature- and 
frequency-independent prefactors, since the RSRG is not 
expected to capture these correctly. In Eq. (|13afl . we see 
that the delta function constraint restricts attention to 
C = 0, i.e., to those pairs of states that are about to be 
decimated at energy ui, leading to Eq. I|13bll . 

Having obtained an expression for [Imif R (s, s'; w)]dis.i 
we now insert it into Eq. (jllafl . We can immediately 
evaluate one of the integrations over the parameters s and 
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s' giving the Laplace transform of P(0, \s — s'\;u>). Up 
to a constant of order unity, the remaining integrations 
J dbds oc A, so that we have 



[fiMJdis. 



cr - 



il tanh(w/4T) 
w log 2 fi /w 



(14) 



where <7o is a dimensionful prefactor that is difficult to 
estimate (in part due to the uncertainty in the prefac- 
tor of rifj); a very rough estimate along the lines of our 
estimate of fio gives ctq ~ e 2 eFAo£/£loa. 




FIG. 2: Plot of [o"i(u;)]dis.(in arbitrary units). Here, to is 



measured in units of 10 
and fi = 10 13 s _1 . 



and we have chosen T = lOmK 



This formula (plotted in Fig. applies at ui <C f2 
and includes only the contribution due to extended im- 
purities. It exhibits a peak in u> which becomes larger 
and narrower and moves to zero as T — > 0. Also, 
[<7i(0)]dis. = 0. This contribution would occur in parallel 
to those arising from pointlike impurities^, and would 
therefore be most noticeable at low temperatures where 
the peak is sharpest. For example, according to the cal- 
culations of Durst and Lee^, the microwave conductiv- 
ity due to pointlike impurities should attain a universal 
value at low temperatures and frequencies, in striking 
contrast to the structure of Eq. (|14|) which would exhibit 
strong variations in this regime. At the typical temper- 
ature scales of the Turner et al experiments^, however, 
the extended impurity contribution would be particularly 
broad and small making a quantitative comparison dif- 
ficult. To isolate the contribution considered here, we 
make two suggestions: (1) This contribution is largest at 
u> — > with T < ui so that the argument of the tanh 
in Eq. I|14|) is large. (2) By deliberately introducting 
extended impurities via ion irradiation^ or by perhaps 
introducing disorder away from the CuC>2 planes, this 
contribution will be enhanced. 



V. CONCLUDING REMARKS 

In an effort to understand the electronic properties of 
disordered d-wave superconductors, we have computed 



the contribution to the microwave conductivity arising 
from extended impurities. This contribution has its ori- 
gin in the hybridization of ZE Andreev bound states oc- 
curring near each such extended impurity. To capture 
this hybridization, we have applied a real space renor- 
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on lattice models to the continuum 
Hamiltonian Ha- From a theoretical point of view, this 
generalization is possible because the underlying super- 
symmetry of Ha^^^ implies that there is an approxi- 
mate zero-energy state near each zero of A(s) for every 
realization of the random pair potential, providing a nat- 
ural basis for studying the low-energy properties of Ha 
via such tunneling corrections. We found that the singu- 
lar density of states associated with these Andreev bound 
states found in Ref. leads to a singular contribution to 
the microwave conductivity and discussed the possiblity 
of observing this experimentally. Finally, we note that 
whereas the low-energy quasiparticles studied here rely 
on the d-wave symmetry of the pair potential, they are 
not associated with the nodes of the pair potential. 
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APPENDIX A: SINGLE-PARTICLE GREEN 
FUNCTION 

In the present section, we apply the RSRG technique 
to the calculation of the Green function of Ha averaged 
over all realizations of the random pair potential. These 
results apply more generally to any system in which Ha 
arisesiSiSLSi. This calculation is similar to the one pre- 
sented in the main body of the text and obtains a result 
which has already been obtained by Balents and Fisher 21 . 
Nevertheless we present it here for the sake of clarity and 
completeness, as well as to highlight the strengths and 
weaknesses of the RSRG method. Indeed, the fact that 



we are able to reproduce (up to numerical prefactors) 
the exact results of Ref. |2lJ lends support for this tech- 
nique. The disorder-averaged Green function G(x~y; iuo) 
is given by 

G(x-y;iu) = [G(x,y;iu))] diB , (Al) 
(iuj- H A )G(x,y;iu) = 5(x-y), (A2) 

where the square brackets in Eq. IjAlfl denote the average 
over all realizations of the pair potential. It is useful to 
express G(x, y; ioS) in terms of the normalized eigenstates 
ip n of H A : 



G(x,y;iuj) = ^ 



1pn(x)lp n (y) 

iuj — E n 



(A3) 



Before proceeding with the RSRG calculation, we em- 
phasize that it does not correctly obtain overall prefac- 
tors associated with the contributions to G(x — y;iu>). 
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Nonetheless, it obtains in a simple way the correct (non- 
trivial) asymptotic behavior of G{x — y; iuS) as well as the 
correct matrix structure. In the following, these are the 
features we are primarily interested in obtaining. 

The RSRG simplifies the computation of G(x — y; ito) 
in two distinct ways. Firstly, by decimating until Q, = to, 
we eliminate all strongly coupled zero-energy (ZE) states. 
Since the remaining ZE states are weakly coupled, it 
is plausible to approximate the exact low-energy eigen- 
states ip n {x) by perturbed ZE states. Secondly, due to 
the fact that P(C, l\ is very broach, certain rare pair- 
potential configurations dominate the average over all 
random pair-potential configurations. Here, we focus on 
two leading contributions; analagous contributions were 
considered in Ref. [l6| in the context of the related prob- 
lem of computing the single particle Green function for a 
one-dimensional random hopping problem. However, we 
obtain slightly different results for the second contribu- 
tion. The first we denote by G\ (x — y; ilu) . This conribu- 
tion is due to pair potential configurations such that the 
coordinates x and y are near sign changes that are about 
to be decimated (owing to the fact that the pair potential 
connecting them is relatively small) in the sense discussed 
in Sec. IIIII The second contribution, G2(x — y; ito), is due 
to pair potential configurations having only one of x or 
y be about to be decimated along with a ZE state at a 
third site. Let us first, however, discuss G\(x — y; ito). 

A A 




x and the + sign for the ZE state ^ y at y: 

* x ex e*- (A4) 

*^Qe # «, (A5) 

where we have suppressed unimportant normalization 
factors and displayed phase factors in these expressions 
which were suppressed in Eq. Q . These factors are cho- 
sen to satisfy expi((/) x — <f) y ) — —i so that t n in Eq. (jSJ) 
is real and positive; one must account for these prop- 
erly to obtain the correct phase of G\(x — y; iuo). Under 
the influence of the tunneling matrix element t, (which is 
small by virtue of the decimation procedure) , these states 
form symmetric (S) and antisymmetric (A) combinations 
ips,A — -y^i^x^^y) with energies ±i(= ±w/2 since they 
are about to be decimated). As in Sec. HVI we take these 
approximate eigenstates to be the lowest energy states 
appearing in the sum Eq. I|A3() (i.e., n = S, A). Thus, 
we find (including the contribution due to configurations 
with A < between x and y) 

Gi(x-y;iuj) cx ct 3 -^ 

oo 

x^t-ir+'A-" 2 * 2 ! 1 ^!^, (A6) 

n=l 

where T = logf^o /u> and by "oc"we specifically mean 
equality up to a real and positive prefactor not captured 
within the RSRG approach. We emphasize that the 03 
structure of this contribution comes from correctly iden- 
tifying the spinor structure of the ZE states in Eq. (|A4|) 
and Eq. i|A5|l . The infinite sum arises from the inverse 
Laplace transform requried to compute P(0, \x — y\;u>). 



FIG. 3: Class of pair potential configurations providing the 
contribution to Gi(x — y\iu). The smallness of A between x 
and y is meant to indicate that the ZE states associated with 
these sites are about to be decimated. 

The rare pair-potential configuration that leads to 
Gi(x — y; iuj) is the same as was considered in Sec. lIVI for 
[lmK R (s; w)]dis. and is depicted schematically in Fig. |3 
the fact that both x and y are about to be decimated at 
energy u> is depicted pictorially by having A be relatively 
low between x and y. Of course, it is not necessary that 
the sites x and y be nearest-neighbor sign changes as de- 
picted in Fig. |21 but merely that they are nearest neigh- 
bors after having decimated until to. The relative num- 
ber of such configurations is given by n w P(0, \x — y\;u)); 
to obtain G\(x — y;icu) we simply multiply this by the 
magnitude of the contribution given by Eq. (|A3p . For 
concreteness, let us assume that A > on average (as in 
Fig. between x and y (y > x); once we have computed 
this it is straightforward to obtain the A < configu- 
rations. The ZE states associated with these points are 
given by Eq. Q with the — sign for the ZE state ^> x at 



A M 




FIG. 4: Second-order class of pair potential configurations 
contributing to Gi(x — y;iuS). The ZE states at sites z and 
y are about to be decimated due to the fact that they are 
strongly coupled (indicated by the relative smallness of A 
between them) but the states at sites x and z may still be 
weakly coupled. 

Next, we turn to the evaluation of G2{x — y; ioS). The 
type of pair-potential configuration we have in mind is 
depicted in Fig. |@J). Again, these states do not have 
to be nearest neighbor sign changes but merely nearest 
neighbors after having decimated until to. The ZE states 
at x and z are weakly coupled (by a matrix element S), 
owing to the relative largeness of A between them. The 
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ZE states at z and y are strongly coupled with matrix 
element t = cj/2 and are thus about to be decimated. 
Let us consider the Hamiltonian H for our system in the 
subspace of these three ZE states: 



H 




denote the amplitude of each ZE state contribution asso- 
ciated with each site. Since the ZE state spinors at x and 
y are each of the form of Eq. (|A5J) . the term appearing 
in the sum has the form (switching to a slightly different 
notation and neglecting normalization factors as usual) 



(A7) 



i)o{x)ij)o{y) 



ts 



1 



The eigenvalues of this matrix are and ±\/t 2 + S 2 ; as 
we are interested in low energies, for the purposes of the 
sum in Eq. (|A3|) we proceed by keeping only the state ipo 
having eigenvalue in Eq. i|A3() . which has the explicit 
form 



1 -t 

y/S 2 + t 2 itu M 1 



(A9) 



where we note that an extra — sign in Eq. I|A9[) came 
from the phase factors in Eq. (|A4|) and Eq. I|A5|I . By 
including the possibility that A < between z and y, it is 
straighforward to see that the sum of these contributions 
(A8) is proportional to the unit matrix. Since it is about to 
be decimated, t = lu/2. The other matrix element S may 
be labeled by £ (as usual) via S — f e~ ? . Thus, we have 
Of course, the elements of the column vector in Eq. (|A8|) [neglecting 6 <C t in the denominator of Eq. (|A9|) ] 



4>o 



\M 2 + 1 2 




G2{x - y; ilu) a^ — 
iui 



dz 



d(P((,x-z;u;)P(0,z-y;Lo), 



(A10) 



r 



where tro is the unit matrix. Taking the Laplace trans- 
form of both sides of Eq. lAlOfl and evaluating the in- 
tegral over C, we find that the Laplace transform of 
G 2 (x - y\ iu) is 



G 2 (y;iuj) oc 



cr 1 



1 



1 L 2 sinhTy/y shinTy 7 ^ + \fy coth Ty/y ' 

(All) 

with T = logOo/^ as before. To obtain G 2 (x — y;iuj), we 
must evaluate the inverse Laplace transform of Eq. (|A1 l|l . 
To do this, we use the standard technique of identifying 
the poles y n of Eq. (|A1 lfl ; these occur when sinh Ty/y^ = 
0, i.e., y n = —n 2 ir 2 /r 2 (we must exclude the case n = 0). 
Multiplying the associated residues by exp(\x — y\y n ) and 
summing over all n yields the result 



Our final result for G(x — y; iuj) — G\{x — y; iui) + G 2 {x — 
y;iui) agrees with the exact results of Ref. (|21r) . as can 
be seen by taking the M — > (M being proportional to 
[A(s)]di s in our notation) limit of Eq. (5.30) in that pa- 
per. In the present context, this form for G(x — y; iuo) 
indicates that the Andreev bound state Green function 
decays exponentially as a function of |a; — y\ but with 
a characteristic length scale log 2 VIq/oj that diverges as 
lj — > 0, indicating the derealization of low-energy quasi- 
particles along a trajectory. 



G 2 (x-y^)oc^Y,^- nVlX - yWr2 

71=1 



(A12) 
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